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Abstract 

We present detailed heat capacity measurements for Ba(Fei_ a; Co a ;)2As2 and Ba(Fei_ 2; Ni, c )2As2 
single crystals in the vicinity of the superconducting transitions. The specific heat jump at the 
superconducting transition temperature (T c ), AC P /T C , changes by a factor ~ 10 across these series. 
The AC P /T C vs. T c data of this work (together with the literature data for Ba(Feo.939Coo.o6i)2As2, 
(Bao.55K .45)Fe2As2, and (Bao.6K .4)Fe2As2) scale well to a single log-log plot over two orders of 
magnitude in AC P /T C and over about an order of magnitude in T c , giving AC P /T C oc T c 2 . 

PACS numbers: 74.25.Bt, 74.62.Dh, 74.70.Dd 
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The discovery of superconductivity in F-doped LaFeAsO^ and K-doped BaFe2As2^ com- 
pounds resulted in a large number of experimental and theoretical studies of the materials 
containing Fe-As layers as a structural unit. The understanding that transition metal (TM) 
substitution for Fe in the (AE)Fe2As2 series (AE = Ba, Sr, Ca) could be used to both sta- 
bilize superconductivity^^ and simplify growth while improving homogeneity, makes the 
(AE)(Fei_ a; TM a ;)2As2 series (and in particular Ba^ei-^TM^As^ 1 ^ 1 ^ 1 ^) a model system 
for studies of physical properties of Fe-As based materials. 

Despite the large experimental and theoretical effort to understand the nature of super- 
conductivity in Fe-As based materials, there are still a number of open issues that nave not 
been well resolved. From the experimental point of view, for number of properties, either the 
spread in the data is large, or systematic sets of data are absent. Temperature-dependent 
specific heat measured through the superconducting transition is one such property, even 
though it is considered to reflect of the superconductivity mechanism. Apart from the 
functional dependence of the specific heat (C P (T)) below the superconducting transition 
temperature (T c ), the study of which often has complications caused by the need of careful 
substraction of the normal state contributions, the size of jump in C p at the supercon- 
ducting transition is known to depend on the details of the superconducting state^^^, 
as judged and modeled by its deviation from the isotropic, weak coupling, BCS value of 
ACp/jTc = 1.43 (7 being the normal state electronic specific heat). 

The Ba(Fei_ x Coa;)2As2 and Ba(Fei_ x Ni IE )2As2 families of materials share very similar 
and complex x — T phase diagrams^ 1 ^: initially on Co(Ni) -doping the critical temperature 
of the structural/antiferromagnetic transition decreases, with a separation in critical 
temperatures between these two transitions,— 1 ^ then, above some critical concentration, 
superconductivity is observed, apparently in the orthorhombic/antiferromagnetic phase. 
At higher Co(Ni) concentrations the structural/magnetic transitions are fully suppressed, 
whereas superconductivity appears to persist in the tetragonal phase up to 1 ~ 0.14 for 
TM = Co and x ~ 0.08 for TM = Ni. In this work we study the evolution of the specific 
heat jump at T c with TM = Co and Ni doping, for different concentration of TM. We 
examine the whole superconducting dome (both in orthorhombic/antiferromagnetic and 
tetragonal low temperature phases of the x — T diagram) to gain insight into the details of 
the superconducting state in these materials. 
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Single crystals of Ba(Fei_ a; Co 3 ;)2As2 and Ba(Fei_ a; Ni a ;) 2 As2 were grown out of self flux 
using conventional high-temperature solution growth techniques.-^ 7 - Detailed description of 
the crystal growth procedure for this series can be found elsewhere.- 1 ^ The samples are plate- 
like with the plates being perpendicular to the crystallographic c-axis. The heat capacity 
data on the samples were measured using a hybrid adiabatic relaxation technique of the 
heat capacity option in a Quantum Design, PPMS instrument. Part of the C P (T) data for 
Ba(Fei_a;Co x )2As2 were presented, but not analyzed in detail, in Ref. [7. 

It has to be noted, that both the superconducting transition temperatures and the 
upper critical fields in these materials are rather high 7 , thus making a reliable estimate 
of the normal state electronic specific heat, 7, difficult, especially bearing in mind that 
for approximately half of the samples in this study superconductivity coexists with an 
antiferromagnetic long range order. For this reason we are limited to the experimental 
determination of AC P /T C , rather than the more traditional quantity AC P I^T C . Due to 
finite widths of the superconducting transitions, AC P /T C and T c values were determined 
from plots of C p /T vs. T using an " isoentropic" construction (Fig. [U(b), inset). So defined 
values of T c may_ 
batches in Refs. 
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3e slightly smaller than those reported for the samples from the same 
lip] ] in which different criterion was used. 



Temperature dependent heat capacity data for Ba(Fei_ a .Co x ) 2 As2 and Ba(Fei_. r Ni, 2 .)2As2, 
plotted as C p /T vs. T 2 are shown in Figs. [Qa),(b). The jumps associated with the super- 
conducting transitions are seen for all concentrations presented. Fig. [T](c) presents all of 
the C P (T) data together showing that the spread of the of the data above the supercon- 
ducting transitions is small, within 5-6%, consistent with simple sample weighting errors. 
In addition, Fig. [Qc) shows that both Co and Ni, added in these small amounts, are small 
perturbations to the BaFe2As2 system and do not significantly change the background C p . 

The values of the specific heat jumps at superconducting transition in the form of AC P /T C 
are plotted as a function of TM = Co, Ni concentration in Fig. El The values of T c as a 
function of x are displayed on the same plots as well. It is remarkable that (i) AC P /T C values 
change by as much as a factor of ~ 10 for the samples within the series; (ii) the shapes of 
the C p /T c vs. x curves for both series appear to be related to the shapes of the respective 
T c vs. x superconducting domes. 

These data can be plotted as AC P /T C vs. T c (Fig. [3]). It is curious, that all the data (both 
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for " underdoped" and "overdoped" parts of the superconducting dome) collapse rather well 
onto a single curve. A data point 8 from another group's work on Ba(Feo.939Coo.o6i)2As2 is 
consistent with our data, moreover, our previous result^ on Bao.55K a45 Fe2As 2 as wel 
the very large C p jump at superconducting transition reported for Bao.6K .4Fe2As2 
follow the same trend (Fig. [3]). All the data in this figure can be fitted by a straight line 
with a slope n 2. 
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There are several possible ways to address such a remarkable behavior. 

Inhomogeneities: it has been known that in a few cases (e.g. K-doping in BaFe2As2 
samples^ 8 - 1 ^ the resulting samples had a distribution of dopant concentration, resulting in 
a broadening of the phase transitions. It appears to be less the case for TM doping: the 
wavelength dispersive x-ray spectroscopy does not show unambiguous, beyond the instru- 
ment error bars, distribution of TM dopant^^, for the more studied Ba(Fei_ x Co x )2As2 
series, the low field susceptibility does not show variation in either transition width or mag- 
netic flux expulsion, 7 Meissner screening is homogeneous,—, structural and antiferromagnetic 
transitions for the intermediate concentrations are reasonably sharp 7 - 1 ^ 1 ^, and, for the phase 
diagrams and several other properties, results from the different groups appear to be very 
close. 7 ' 8 This being said, in an oversimplified model, inhomogeneities can be modeled by a 
uniform distribution of superconducting transition temperatures within some temperature 
range. For an aggressive, ±10% of T c spread, the width of the distribution (using x = 0.074 
Co data as a starting point), the apparent AC P /T C jump will be approximately factor of 
3 (but not ~ 10) smaller than the initial one, thus suggesting that inhomogeneity is not 
the sole reason for the observed behavior. In addition, the fact that the data in Fig. [3] are 
linear, with a slope n ~ 2, over two orders of magnitude in AC P /T C and over about an order 
of magnitude in T c on a log-log plot argues against an artifact caused by an uncontrolled 
spread in composition and may imply some more profound physical mechanism. 

Significant change of the density of states within small, ~ 10% range TM = Co, Ni dop- 
ing, with j(x) or DOS(x) having dome-like shape centered at x values corresponding to the 
observed maximum in T c : in this case AC p j^T c could be close to constant or change insignif- 
icantly, whereas strongly x-dependent AC P /T C will be observed. Although, as mentioned 
above, reliable experimental data on normal state electronic specific heat as a function of x 
are not available, band structure calculations on pure BaFe2As2r^^ 7 - do not suggest such 
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a significant and non-monotonic change of the density of states for small TM = Co, Ni 
concentrations (a sharp local maximum in 7 would be required at optimal doping). 

More physical reasons to observe significant change (decrease in comparison with 
isotropic, weak coupling, BCS case) in the heat capacity jump at T c could be associated 
with the effect of paramagnetic impurities,— multi-band superconductivity,— ^ and/or 
effects of a normal state pseudogap.— Although each of these possibilities is plausible and 
exciting, it seems hard to construct a simple picture that will accommodate the observed 
dome-like, almost symmetric, AC P /T C vs. x behavior within a single one of this models 
(unless a more complex case, in which left and right, "underdoped" and "overdoped", parts 
of the AC P /T C vs. x dome, are explained separately by different mechanism, is considered). 

To summarize, approximate scaling of AC P /T C with T c was observed for 
Ba(Fei_ x Co x ) 2 As 2 and Ba(Fei_ x Ni x ) 2 As 2 single crystals. The reason for such scaling is 
not clear in this moment: if extrinsic (chemical inhomogeneities) it cannot be ignored in 
interpretation of other, detailed, experiments on these materials, if intrinsic, more work is 
required to elucidate the reason for this apparent scaling. 
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FIG. 1: (Color online) Temperature-dependent specific heat of (a) Ba(Fei_ x Co x .)2As2, (b) 
Ba(Fei_ x Nia;)2As2 single crystals plotted as C p /T vs. T. Inset to panel (b): enlarged C p /T vs. T 
plot near the superconducting transition for Ba(Feo.954Nio.o46) lines show how T c and AC P /T C are 
estimated. Data in panels (a) and (b) are shifted by a multiple of 50 mJ/mol K 2 along the y - axis 
for clarity. Panel (c): data for both series plotted as C p /T vs. T 2 without shifts. 
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FIG. 2: (Color online) AC P /T C (circles, left axis) and T c (asterisks, right axis) as a function of 
concentration, x, Ba(Fei_ a: Co a ;)2As2 (upper panel) and Ba(Fei_ a; Ni a ;)2As2 (lower panel). 
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FIG. 3: (Color online) AC p /T c vs. T c for Ba(Fei_ a; Co a ;)2As2 (circles) and Ba(Fei_ a; Ni a ;) 
(stars) plotted together with literature data for Ba(Fe .93gCoo.o6i)2As2 [8j, Bao.55K .45Fe2As2 
and Ba .6Ko.4Fe 2 As2 QQ- 

Dashed line has a slope n = 2 and is a guide for the eye. 
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